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A confined system of non-interacting electrons, subject to tfie combined effect of a time-dependent 
potential and different external chemical-potentials, is considered. The current flowing through 
such a system is obtained for arbitrary strengths of the modulating potential, using the adiabatic 
approximation in an iterative manner. A new formula is derived for the charge pumped through an 
un-biased system (all external chemical potentials are kept at the same value); It reproduces the 
Brouwer formula for a two-terminal nanostructure. The formalism presented yields the effect of the 
chemical potential bias on the pumped charge on one hand, and the modification of the Landauer 
formula (which gives the current in response to a constant chemical-potential difference) brought 
about by the modulating potential on the other. Corrections to the adiabatic approximation are 
derived and discussed. 
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I. INTRODUCTION AND SUMMARY 

The flow of a dc current in response to a slowly- varying 
time-dependent potential operating on an un^biased sys- 
tem is termed "adiabatic charge pumping" .EJ This phe- 
nomenon, first considexeiL in Ref. ^, has at±ija|C±ed re- 
cently much thcorcticalJLj and experimentallljilZl inter- 
est. In general, "adiabatic pumping" occurs when the 
charge transferred across a boundary during a single pe- 
riod of a certain modulating potential is independent of 
the modulation frequency. This process is adiabatic in 
the sense that the periodic potential varies very slowly 
in time, such that its frequency uj is smaller than any 
characteristic energy scale of the_filectrons. 

Eighteen years ago, ThoulessEJ has shown, using the 
adiabatic approximation, that the ground state of an infi- 
nite one-dimensional (ID) system of non-interacting elec- 
trons subject to a slowly moving periodic potential can 
support a dc current. Later theoretical investigations of 
quantum pumping in confined nanostoictures have em- 
ployed the result derived by Brouwer J3 which gives the 
pumped charge in terms of the time-dependent scattering 
matrix related to the modulating potential. The deriva- 
tion presented by Brouwer is based on the analysis of Ref. 
p^ , which, in turn, utilizes an expansion in the ampli- 
tude of the modulating potential in the context of time- 
dependent scattering theory. Nevertheless, it is_generally 
accepted that the formalism of Biittiker et al&B and the 
resulting Brouwer formula are validcj even for large am- 
plitudes, as long as the adiabatic approximation holds. 

This paper is devoted, among other issues, to the ex- 
ploration of this point. We consider a spatially-confined 
system of non-interacting electrons, connected by leads 
(denoted by a) to electronic reservoirs which are kept 
at various chemical potentials, fia- The system is also 



subject to a slowly-varying periodic potential. We de- 
rive an expression for the instantaneous current flowing 
in this system, allowing for arbitrary strengths of the 
modulating potential, using time-dependent scattering 
states. The formalism is based on an iterative solution 
of those states, in the adiabatic approximation, in which 
the temporal derivative of the scattering potential (and 
the scattering states) is the small parameter. The formal 
derivation is summarized in Sec. II, where we obtain the 
instantaneous current in the lowest-order adiabatic ap- 
proximation. This current is averaged over a single pe- 
riod of the modulating potential. In this way we obtain 
the effect of the chemical potential bias on the pumped 
charge on one hand, and the modification of the Landauer 
formula caused by the modulating potential on the other. 
We find a new expression for the charge pumped through 
an un-biased system [see Eq. (|^) below], which is par- 
ticularly useful in cases where the modulating potential 
operates on the entire nanostructure, Ej and is spatially- 
dependent. The Appendix includes the derivation of the 
next order correction to the current. 

We investigate the lowest-order expression for the cur- 
rent in Sec. Ill, confining ourselves for simplicity to a sys- 
tem connected to two reservoirs. Our results there can 
be summarized as follow. In a two-terminal structure, 
the current averaged over a single period, r, consists of 
two parts. The first, denoted by /pump, fiows even when 
the system is un-biased, (but is modified by the presence 
of the chemical potential difference). To the lowest order 
in the adiabatic approximation, that current reads 
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where £ and r denote the left and the right leads, respec- 
tively, and fe,r{E) = (e<-^-t^i.r-)/kBT _,_ ^re the Fermi 
distributions in the reservoirs connected to the left and 
right leads. In Eq. (|]), |/3m) = x^m ^^'^ instantaneous 
scattering state at time t, excited by an incoming wave 
in channel m of lead /3, and V is the temporal derivative 
of the time-dependent scattering potential. Interestingly 
enough, the matrix element appearing in Eq. can 
be written in terms of the instantaneous scattering ma- 
trix, reproducing the Brouwer formula,^ and vindicating 
its use for arbitrary amplitudes. The second part of the 
current, denoted by /bias, flows only when the system is 
biased, i.e., ^ ^r- When this current is integrated over 
a single period, r, of the modulating potential, one ob- 
tains the Landauer formula, modified by the modulated 
potential in two ways: (i) The transmission coefficient 
is the instantaneous one, T*, averaged over a single pe- 
riod of the time-dependent potential; (ii) There appears 
a correction to the Landauer expression, which is related 
to the temporal derivative of the modulating potential 
and the ensuing instantaneous scattering matrix. In the 
simplest case in which the two leads are single-channel 
ones, /bias takes a particularly simple form. 



/bias — 



dt 



dE 
2 dt dE 
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in which -0* is the instantaneous Friedel phase (i.e., the 
transmission phase) of the nanostructure. The first term 
in Eq. (^ yields the Landauer formula for the present 
case. The second term there is a correction, which is 
discussed in Sec. Ill, using a simple example. 



II. TIME-DEPENDENT SCATTERING THEORY 
IN THE ADIABATIC APPROXIMATION 



scattering solutions can be defined un-ambiguously. This 
confined region is attached to leads, numbered by the in- 
dex a, and each lead is connected to a reservoir having 
the chemical potential fXa- The Hamiltonian Tio con- 
sists of the kinetic energy. We use the adiabatic approx- 
imation, requiring that the characteristic inverse time- 
constant, l/r, which describes the time dependence of 
V, is smaller than any characteristic energy scale of the 
electrons. For a simple oscillatory potential, r ~ 211/10. 

As in the usual scattering treatment, we denote the 
incoming wave with energy E in lead a by ■u;~„ , where n 
is the transverse mode number. This wave is a solution 
of the free Hamiltonian, 



ino{r)~E)w-^^{r)^0, 



(4) 



and is normalized such that it carries a unit flux. The 
scattering solution of the full Hamiltonian, excited by 
this incoming wave, can be written in the form 



Xanir, t) ^ W^J^X) + Xa«(r, t). 



(5) 



The time dependence of the scattered wave function, 
XQn(r, t), is expected to have the same characteristic time 
scale as V . For example, when the modulating potential 
is oscillating in time, x contains all harmonics of the fre- 
quency LO. 

The scattering solution '^an should satisfy the time- 
dependent Schrodinger equation. 



.5*a„(r,i) 

« = W r, t)^'a„(r, t). 

at 

Inserting (||) into (^), using Eq. (Q), we find 

.9xan(r,i) 



(g\E)^ Xcnir, t) = Vir, t)w-Jr) - 



dt 



(6) 
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In the first part of this section we solve for the time- 
dependent scattering states iteratively, using the adia- 
batic approximation. We then use those scattering states 
in the second part, to obtain the current. The formalism 
presented below borrows from the derivations in Refs. 

and extended to include the effect of a time- 
dependent scattering potential. 



A. Time-dependent scattering states 

We consider a ballistic nanostructure of arbitrary ge- 
ometry, which consists of a nanostructure connected to 
several electronic reservoirs. This system is described by 
the Hamiltonian 



nir,t)^noir) + Vir,t), 



(3) 



where the scattering potential V{r,t) is assumed to be 
confined in space, so that asymptotic behaviors of the 



where G* is the instantaneous Green function of the full 
Hamiltonian. such that 



E -n{r,t))G\E;r,r') ^ Sir' - r). 



(8) 



We now solve Eq. (^, using the adiabatic approxi- 
mation: the temporal derivative appearing on the right- 
hand-side of that equation is regarded as a small correc- 
tion (of order 1 /t) , and the equation is solved iteratively. 
The zeroth-order is just the instantaneous scattering so- 
lution, which we denote by Xan, 

xLir) ^ w-nir) + j dv'G\E-v,v')V{v',t)w^^{v'). (9) 

The instantaneous scattering state x*an is the solution of 
the instantaneous Schrodinger equation, with energy /?, 



i?-H(r,i))xL(r) 



0. 



(10) 
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Turning back to Eq. one finds that to first order in 
the time-derivative the scattering solution reads 



in which is the time-derivative of the instantaneous 
scattering state. Hence, to first order in the adiabatic 
approximation the time-dependent scattering states are 
given entirely in terms of the instantaneous solutions 
(namely, Xam G'*(£'; r, r')) of the problem at hand. 
One notes that the adiabatic solutiop-,(|ll|) of the scat- 
tering state is analogous to Thoulesala solution for the 
ground-state wave function in his model. In the Ap- 
pendix, we discuss the corrections to the lowest order 
adiabatic approximation. 



B. The current 

Here we outline the derivation of the current in the 
scattering states formalism, as developed, e.g., in Refs. 
p2| , p3| , and |2j. One writes the field operator of the elec- 
tron, 5'(r,t), in terms of the scattering states as 



-lEt 



(12) 



in which dan destroys an electron incoming in channel n 
of lead a. The thermal average of the latter operators is 
given by the Fermi distributions of the various reservoirs, 
such that 

{ai„iE)aa'n'iE')) = 27rSiE - E')San,a'n' faiE), (13) 

where fa{E) is the Fermi distribution in the reservoir 
connected to the a lead. With these definitions, the ther- 
mal average of the current density operator becomes 



5:/.(£;)xL(r,t)%M, (14) 



dv 



where e stands for the negative electron charge. 

It is convenient to evaluate this quantity when r ap- 
proaches oo in lead (3 (which will be denoted by r cxd/3), 
and then to integrate the current density over the cross- 
section of that lead (noting that the incoming and out- 
going waves are normalized to carry a unit flux). In so 
doing, we may take advantage of the asymptotic prop- 
erties of the instantaneous quantities x*an s-i^d G* {E) , as 
documented in Refs. |23 and 



G'{E-v,v% 



-♦00/3 



xL(l")|r^oo/3 = (5a/3W„„(r) + 



(r)5, 



f3m,an • 



(15) 



Here, is the outgoing wave in channel m of lead /3 
and S*^^ is the matrix element of the instantaneous 



scattering matrix. As a result, the current flowing into 
lead P is given by 



Xa„(r,t) = xL(r) - ^ J dr'G'{E;r,r')xL{r'). (11) ^t) = ^ J 



dE ^ 



f37n,an I 



25ft (, 



with 



Upm.an = / drx^„ (r)xL (r) ■ 



(16) 



(17) 



The result (|T^) for the time-dependent current entering 
into lead (3 of the nanostructure holds for a general biased 
system, whose various terminals have different chemical 
potentials (as long as the time-dependence of the peri- 
odic potential is slow enough). It is therefore interesting 
to consider charge conservation, using that result. In- 
deed, summing (Qq) over all leads, we obtain 



E 





(18) 



which shows that when the total current entering the sys- 
tem, X]/3^^(0; is integrated over a single period of the 
modulating potential, the result is zero, i.e., the charge 
per period is conserved. In deriving the result ( p^ we 
have employed (i) the unitarity of the instantaneous scat- 
tering matrix, Xi/J™ Sfm.anSpm.a'n' = ^a«,a'«'p|nd (11) 

the following property of the scattering matrixOEJ 



/3m 



(19) 



Equation (|16|) can be considered as a generalization of 
the Landauer formula, extended to include the effect of 
a time-dependent potential, in the adiabatic approxima- 
tion. The new ingredient is the quantity Uan,i3m, Eq. 
(p!?!). This quantity can be expressed in terms of the 
temporal derivative of the scattering potential. 



/3m, an 



rfr(-%^)T>(r,t)xL(r). 



(20) 



To prove this, we take the temporal derivative of Eq. 
(|l0|), and use Eq. (||), to obtain 



(r) 



j dr'G\E;r,r')V{r',t)xUr'). 



(21) 



We insert this expression into Eq. (p7[), and carry out 
one of the spatial integrations using 



J dr'G\E;r,r')xLir') = - 



dE ' 



(22) 



which follows directly by differentiating Eq. ([T^) with 
respect to the energy, and using Eq. (^) and the sym- 
metry of the Green function G*(£';r,r') = G*(i?; r',r). 
This produces the result (20). 
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III. THE TWO-TERMINAL SYSTEM 

Let us now confine ourselves to a nanostructure con- 
nected to two terminals, with left {£) and right (r) leads. 
Then we can use Eq. to write the current entering 
the system from the left terminal in the form 



dE 



+ [h{E) + U[E) 



h{E)- U{E))(\Sl 



of TT* 

£m,rn £m,rn 



(23) 



An analogous expression holds for Ir{t). The net cur- 
rent flowing in the system during a single period of the 
modulating potential then consists of two parts, 



pump 1 



(24) 



where the first, /bias) flows only when the system is bi- 
ased, whereas the second, /pump, is established by the 
time-dependent potential (though it is affected by the 
chemical potential difference, when the latter is applied). 
Using Eqs. ( ^ and (po|), the pumped current, /pump, 
takes the form 



/pump ^ 



dE 



-(h{E)+fr{E) 

n V 



(25) 



For the biased current we find 



dt f dE f 



/bias ^ e (t — / — 



E 



2|'5'rm.^ril 



In 



u* 



s„ 



-[ME)^fr{E) 
TT V 



U* 



rm.rn^ rm.rn 



(26) 



The pumped part of the current, /pump, can be written 
in terms of the temporal derivatives of the instantaneous 
scattering matrix. In order to vindicate this statement, 
we start from the asymptotic form for x*am ^IQ- (jH), for 
the wave going from lead a into lead (3. Noting that 



,(r ^ 00, /3) = w^Jap + (1 - 5ap) ^ w+^^, (27) 

m 

we conclude that 



5^™,a„ = 1 - <5a,3 - j / drxpm{Y)V{v,t)w^^{v). (28) 



Differentiating this expression with respect to time and 
using Eqs. (y) and (21), yields 



''^ f37n,an 



*X^m(i-)^(i-,i)xL(r)- 



(29) 



Note that this expression is not a matrix element. To 
turn it into an expression involving matrix elements, we 
use Eq. (|l|), by which 



iSlm^c.n = J c^i-y(r,t)x^„(r) 



^ 5^ .anXp'm'i^) 



(30) 



Then, multiplying by S^^^ ^„ and summing over a and n 
yields 



One notes that this identity makes the expression fop; 
/pump, Eq. (|2^) above, to be identical with the BrouwerEl 
formula, in the case where the system is un-biased. 

In the simplest case where each of the leads is a single- 
channel one, the expression for the current takes a partic- 
ularly simple form. In this situation, the instantaneous 
scattering matrix becomes a 2 x 2 matrix, which can be 
parametrized (in the absence of a magnetic field) as 



s- 



(32) 



Here, T* and /?* are the instantaneous transmission and 
refiection, respectively. The reflection phase a* describes 
the asymmetry of the nanostructure. A flnite time- 
dependent reflection phase is a necessary ingredient to 
obtain the pumped current. Finally, if)* is the transmis- 
sion (Friedel) phase. With the parametrization Eq. (|32" 
one finds, using Eqs. (ES) and (pT 



ME)~f,.{E) 
d 



2r 



4(ji(<-*--')f)-i('«'--')i;) 



+ [ME)+friE)) 



_d_ 

dE 



, da' 



(33) 



The charge passing through the nanostructure during a 
single period of the potential, Q, is then 



Q — Qbias Q 



pump- 



(34) 



Here, 



4 



Qbias = — f dt I dE 



d 



(^f,{E)-U{E)) 



2T* 



(35) 



It is seen that the first term here is just the Landauer for- 
mula, with the transmission coefficient averaged over the 
temporal period. The second term forms a correction to 
this result, brought about by the modulating potential. 
The pumped charge is given by 



e 
d 



dt / dE 



x^(/.(i^) + /.(i^))i^*^, (36) 



(see Ref. O), and it vanishes unless a* is time-dependent. 



A. Example— the single-level quantum dot 



Mk = 



2i J sin ka 



Ek-eo + e'''<'{jf + J^)/J' 



(39) 



Let us consider first the pumped charge. At zero tem- 
perature, it is given by 



Q 



pump - 2^ / dt 



da^ 
~dt 
sin ka 



27r 7 "^Vfe - eo + e'^^nXf + X,)|2 
X (X^ — Xr)(eo ^ Ek) + EkQ^r^e — X^X^ 



in which energies are measure in units of J, and 



(40) 



(41) 



Clearly, a comparison between Qp 



and the sec- 



ond part of Qbias is called for. Consider for simplicity 
zero temperature. Then, Qpump is given by the val- 
ues of R^da^/dt a,t E — E-p ± <5/i/e, where J/i is the 
chemical potential difference, and E'f denotes the Fermi 
level. On the other hand, the second term in Qbias 
is given by - J dE[-{S^i/e){d^ f /dE^)T\dil;^ /dt)], and 
hence should be much smaller than the Landauer contri- 
bution, which is proportional to dfi/e. Nevertheless, it 
may be of interest to explore this term experimentally, 
as it is related to the Friedel phase of the nanostructure. 

To further explore this point, we consider the follow- 
ing simple example: a quantum dot, witii-|a single local- 
ized level, coupled to two ideal ID leads.c3 Adopting the 
tight-binding description, we model the two leads con- 
necting the quantum dot to the electronic reservoirs by 
ID chains of sites, whose on-site energies are assumed to 
vanish, and whose nearest-neighbor transfer amplitudes 
are denoted by — J. Thus the energy of an electron of 
wave vector k moving on such a chain is 



Eh = — 2Jcos fco 



(37) 



where a is the lattice constant. The localized level, of en- 
ergy eo, is attached to the left-hand-side lead with matrix 
element —Je, and to the right-hand-side lead with matrix 
element —Jr- The latter two quantities are assumed to 
vary slowly in time, in a periodic way. Our formalism re- 
quires just the knowledge of the instantaneous scattering 
matrix of the system. For the case at hand. 



with 



1+ (^1 Mk 



JtJr 



-1 



Mk 



(38) 



are the time-dependent parameters of the system. Note 
that these two parameters can be thought of as the 'con- 
tact conductances' of the quantum dot. Now imagine 
those to vary in time as follows: Initially, both are equal 
to Xi. Then X^ is increased linearly in time until it 
reaches the value X2, while X^ is being held fixed at the 
value Xi. From that point, X^ is held fixed, while X^ 
increases linearly to the value X2, and so on, see Fig. |l|. 




Xi 

_j 



X2 

_j 



FIG. 1. Tlie periodic temporal evolution of the parameters 
and X,. in tlie parameter plane. 
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It is quite straightforward to find Qpump 
cycle. One obtains 



for such a 



Q 



X2 



pump--/ dX F(X;Xi)-F(X;X2) , (42) 



with 



F(X;Z) = 



sinA;a(eo — Ek + i?feZ) 



(43) 



The resulting charge differs significantly from zero, and 
approaches unity (in units of e) as long as the line of max- 
imal transmission in the Xf-Xr-filane is well within the 
closed orbit forming the period.cl In the present model, 
that line is given by 



X/^ + Xr 



, Ek — ep 

Ek 



(44) 



We will elaborate on this point, which is exemplified in 
Fig. H, in a future publication. For the parameters used 
to produce the results shown in Fig. ^ for Xi = the 
maximal transmission line is contained within the pump- 
ing contour (see Fig. Q), while for its higher values the 
pumping contour shifts away from the maximal transmis- 
sion line. 



. 8 



. 6 



0.4 



0.2 




Xi = 



Xi=l 



Xi=2 



5 10 15 20 



FIG. 2. The pumped charge (in units of e) as function of 
X2, for several values of Xi, indicated on the figure. Here 
ka — 77 — 0.1, and eo ~ —3. 

Let us next consider the time-average that appears in 
the second term of the biased current, i.e., § T^dip* /dt. 
In our example, this quantity becomes 



idiT^'-t- 



dtT'^ , 
J dt 

4sin3 ka{Ek - eo)X^.X^(X£ -I- X^) 



dt- 



|^fc-eo + e''^-(X£-f X,)|^ 



(45) 



On a symmetric periodical curve like that presented in 
Fig. 1^ this integral vanishes. Hence, in such a symmetric 



configuration, there will be no deviation from the Lan- 
dauer formula due to the temporal variation of the Friedel 
phase. One may conclude from this example that when 
the pumping orbit in the plane of the time-dependent 
parameters is not so simple as the one considered above, 
there might be a small correction, related to the asym- 
metry of the orbit. 



IV. CONCLUDING REMARKS 

The motivation of this work is mainly to examine the 
validity of the adiabatic approximation in calculating the 
charge pumped quantum-mechanically through a con- 
fined system of non-interacting electrons, and to study 
the effects of a constant bias. To this end, we have 
derived the time-dependent current generated when the 
system is subject to a slowly-changing modulating po- 
tential, in addition to being connected to reservoirs of 
different chemical potentials. We have found that the 
widely-used Brouwer formula which gives the pumped 
charge in terms of the instantaneous scattering matrix is 
valid for arbitrary amplitudes of the modulating poten- 
tial, as long as the lowest-order adiabatic approximation 
can be employed. 

In the process, we have obtained an alternative ex- 
pression for the pumped charge, Eq. (0), which gives 
it in terms of matrix elements of the temporal deriva- 
tive of the potential between the instantaneous scatter- 
ing states. We have also derived and analyzed the effects 
of the modulating potential on the biased current, and 
showed them to be quite small. Our formalism allows for 
the possibility to obtain systematically the corrections to 
the lowest-order adiabatic approximation. 
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APPENDIX A: CORRECTIONS TO THE 
ADIABATIC APPROXIMATION 



The results presented above have been obtained in the 
adiabatic approximation, in which only first derivatives 
with respect to time are kept. As stated above, this 
means that in the iterative solution of Eq. (|^) only the 
first iteration has been maintained. Here we discuss the 
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contribution of the next iteration. For the sake of simphc- 
ity, we carry out this calculation for an un-biased system 
with two single-channel leads. 

To second order in the temporal derivative, the scat- 
tering solution can be presented in the form 

X(r,i) = x*(r) 



with 



AxL(r') 



dr'G' {E- r, r') (xL (r') + AxL (r')) , (Al) 

-^|dr"|(G*(i^;;r^r")xL(r")). (A2) 

It follows that the quantity C/, Eq. (|l7|), is now modified 
into 

Upa = Upa + AC//3Q,, 

Upo. = / drx^(r)xl(r), 



(A3) 



AC/^„ = / drxMr)AxL(r). 



As a result, we find that the current entering lead /3 con- 
sists of two parts, the leading order in the adiabatic ap- 
proximation, Ip, which has been discussed above, and a 
correction, A/^. Explicitly, 



ip{t) = ip{t) + /^ip{t), 

a 

X '5^p-\SlS + 'i'^{sl*^Upo) 



Alpit) 



'\Upc 



23? 



(A4) 



Let us first verify that the correction /S.Ip obeys charge 
conservation over the entire period. To this end, we sum 
A//3 over j3. Using Eq. (US), we have 



Y,\^Up^\''-2^(sl*^AUp^] 

jdv I dr'[^xMr)£(r)xr(r')xL*(r') 
J J p 

23(xL*(r)|G*(i?;r,r')xL(r') 



(A5) 



We can now employ the fact that the current should be 
conserved upon integrating over the entire period, that 
is, when expression (A5) is inserted into § dt. Then we 
may integrate the-second term by parts. Consequently, 
using the relationE3 



E^aWXa(r')--2QG*(£;;r,r') 



(A6) 



the two terms in (A5) exactly cancel one another. 

Let us now turn to the expression for the correction 
A/^, Eq. (A4), and insert there the condition that the 
system is un-biased, i.e., /g {E ) = f{E). Then, using 
Eqs. © , (|l|), (H), and we find 



Pa\ 



dE 

t* /r.'^ 



V{r,t) 

{G\E;r,r'))nr',t){^^), (A7) 



and 



25ftE5^;Af/0„ = 23 fdr f dr'[[-^)v{r,t) 

a J J 

(dG\E- r,r')^ 



dE 



V{v',t)x'p{v') 



'd'x'pir) 



'd^x'pi^) 



(A8) 



To obtain the last equality, we have made use of the re- 
lations 



dE^ 



= -2 



/*'(^)G*(i?;r',r), (A9) 



and 



(?*(£;; r,r') 
= / dv^G\E-v,v^)V{vi,t)G\E-v^y). 



(AlO) 



Both relations are obtained by taking derivatives of Eqs. 
and (|l^) with respect to the energy and the time. 
Collecting all these terms, we obtain 

Mp{t) = ^ y dEf{E) 

xS^(x^|]/(t)G*(i?) + \v{t)G\E)\x'p). (AH) 

It is satisfactory to note that, again, (upon integrating 
by parts with respect to the energy), the energy integral 
includes the derivative of the Fermi function. Hence, up 
to the second-order in the adiabatic approximation, we 
find 

+ 'i{{x'e\2V{t)G'{E) + V{t)G'{E)\x's))\ (A12) 

To estimate the relative magnitude of the correction (the 
second term above) compared to the leading order one 
(the first term there), consider the part including the 
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second derivative of the potential. Using Eq. (22), this 
part becomes 



dE 



(A13) 



It follows that the correction will be smaller than the 
leading order term by a factor proportional to l/r, due 
to the extra temporal derivative, times the energy deriva- 
tive of the instantaneous scattering state. The latter will 
include the energy-derivatives of the instantaneous trans- 
mission and reflection amplitudes (which appear in x*), 
and possibly a term proportional to the ID density of 
states, i.e., the velocity, (coming, e.g., from the factors 
gifca -^yj^j,;}^ appear in x* of our simple example discussed 
in Sec. III). We may conclude that as long as the en- 
ergy derivative of the instantaneous scattering solution is 
small on the scale of w oc 1/t, the leading order adiabatic 
approximation sufflces to obtain the pumped charge. 
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